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ABSTRACT 



Amplitude and phase relationships between sinusoidal pressure dis- 
turbances in the exit plane of a conical nozzle and mass flux response 
were experimentally determined for the subsonic flow of air. in the 
nozzle. The experimental results were compared with values predicted 
by using the modified asymptotic WKB method. 

Experimental measurements were made with a constant temperature 
hot wire anemometer and a condenser microphone. Disturbance frequen- 
cies from 20 to 250 Hz were investigated for Mach numbers up to .34. 

The experimental values of amplitude factor agree quite well with 
analytical predictions for all Mach numbers tested. Experimental 
values for phase angle differed from those predicted but exhibited a 
similar trend. 
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NOMENCLATURE 



speed of sound 

speed of sound at stagnation conditions 
flow area 

flow area where M = 1 
arbitrary constants 

substantial derivative 

defined in Eq. (24) 
defined in Eq. (29) 
enthalpy 

stagnation enthalpy 



Mach number 
defined in Eq. (62) 
amplitude factor 
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ratio of specific heats (= 1.4 for air) 
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defined in Eq. (25) 
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dimensionless time a t/R 
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dimensionless enthalpy perturbation H / a Q 
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I. INTRODUCTION 



The engineering application of nozzles in the control of fluid flow 
is extremely widespread. It ranges in scope from small fluidic devices 
to large rocket thrusters. In the steady state there exists, for a 
given absolute pressure ratio across the nozzle, a corresponding mass 
flow rate. The introduction of pressure fluctuations at the nozzle 
boundaries affects the flow and certain dynamic characteristics are ex- 
hibited for a given set of conditions, i.e., nozzle geometry, type of 
fluid medium, and the nature of the pressure fluctuations. Because 
Such fluctuations are quite common, especially in systems possessing 
cyclic operating components, a knowledge of the dynamic characteristics 
becomes important. In the past, a quasi-steady relationship between 
mass flew rate and pressure fluctuations was assumed. Thus, an unsteady 
situation was reduced to a sequence of instantaneous steady state cor- 
respondences. The validity of this assumption is questionable. 

Oppenheim and Chilton [3] conducted a literature survey on pulsating 
flow measurement. Most of the investigations noted had to do with the 
effect of pressure fluctuations on the time average of flow. Elrod [1], 
working in connection with a stability analysis of externally pressurized 
thrust bearings, performed a frequency response analysis to predict the 
amplitude and phase relations between the instantaneous flow rate and a 
sinusoidally fluctuating exit pressure in an infinite conical nozzle. 
Although the chosen geometry was extremely simple, Elrod [1] stated 
"It is believed that the essential characteristics of a compressible 
fluid expanding and accelerating from a stagnation reservoir are well 
enough realized for the results to have semi -quantitative practical 
significance.." He further assumed that the flow was adiabatic. 
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frictionless, and irrotational . The analysis using the WKB method is 
asymptotically valid for high frequency number W=coR/a and low Mach 
number. Unfortunately, the numerical results were wrong because of 
a wrong choice of sign for the divergence term of the developing equa- 
tions. 

Chiang, et a!., [2], realizing the generality of Elrod's analysis 
(it can be applied to either converging or diverging flow), corrected 
the choice of sign and extended the region of validity to lower fre- 
quency and higher Mach number by modifying the WKB method. They also 
investigated the problem of pulsating flow in a finite conical nozzle. 

The objective of this thesis was to experimentally investigate the 
mass flux - pressure relationship for pulsating subsonic flow in a 
conical nozzle and to compare experimental results with values pre- 
dicted using the modified WKB solution. 



Elrod [1], in his original analysis, assumed an adiabatic, friction 
free (i.e., isentropic) flow originating in a stagnation reservoir. 
Under these conditions, the flow remains irrotational in accordance 
with Kelvin's theorem and the following governing equations pertain: 



II. THEORETICAL ANALYSIS 



Momentum 



+-C v-p = o 

t>* ? 



0 ) 



Continuity 




( 2 ) 
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which can also be written as 



^ + v-f = o (3) 

Gibb's Equation for a simple compressible substance 



~Tds=dh — L djp=0 (isentropic) (4) 

e r 

For irrotational flow: 

= v-f- (si 

(see Robertson [4], p. 46.) 

Using Eq. (5), Eq. (1) can be rewritten as 

ff z-j-Vf- v-f- (6) 



Stagnation enthalpy is defined as 



H-A+-&- (7) 

z 

Equations (4), (6), and (7) combine to yield 

■f§-= - VH ») 

The scalar product of q with Eq. (1), together with Eq. (4), gives 











DA + aA 
. \)t ** 



(9) 
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which can also be written as 



D H _ J_ a^> 

D* ~ f d± 

by considering a time derivative form of Eq. (4). 
The acoustic velocity is defined by 




Therefore, by the chain rule, Eq. (10) becomes 

J_ JLtL - -L JLg. - aXnP 
CL X ot ~ f d£ d± 



( 10 ) 



( 11 ) 



(12) 



Taking the total derivative of Eq. (12) gives 



D_ T J_ DH " 
f)£ L a. x d* 



D 3>Xn p 

r\ *f* 4- 

V A. 



(13) 



but 



_D_ IAilE - 



p ^ _ 3j 

T±" a* v 0 w 2 



3 

3* 3* 










(14) 



Substituting Eqs. (3), (8), and (14) 'into Eq. (13) above produces 
Elrod's [1] "propagation equation" for stagnation enthalpy: 



D_f± j>H_ 

Di a ‘ D* J 

L eJ 



V J H+ VH- VjUo 



(15) 



For spherically symmetric flow in a conical passage with slight 
flow deviations which are also spherically symmetric, the flow is des- 
cribed by only one spatial coordinate - namely the radius r originating 
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at the apex of the cone. 



Chiang, et al., [2] used linear perturbation theory to express a 
quantity as the sum of a steady-state value and a perturbation. 
Stagnation enthalpy was written as 



H = H+H 



(16) 



where the prime denotes a perturbation term. 

Since the stagnation enthalpy is uniform throughout the mean flow, 
it is apparent that 

VH- 0 07) 

and 



DH 



~ O 



08 ) 



Hence, Eq. (15) may be written, to the first order in small quan- 
tities, as 



_D_ 

D* 



J_ _D j£ 
d x t>± 



= V^H'+VH 7 * Vs&nf 



(19) 



As mentioned previously, spherically symmetrical flow is one- 
dimensional. Chiang, et al., [2] pointed out that for this case 




JL J2_ + ai_ 

D* 7) t ^ dr 



The upper sign pertains to diverging flow and the lower sign is 
related to converging flow. This choice of signs is based upon. the 
fact that for diverging flow the velocity and increasing radius are 



12 



identical in direction. For converging flow, the directions are opposed. 
Therefore, the lower (negative) sign applies to nozzle flows. 

For convenience, Elrod [1] introduced the following dimensionless 
variables: 







(21) 



where a Q is the acoustic velocity at stagnation conditions and R is 
the fictitious radius at which M = 1 for conical flow. In terms of 
these new variables, the "propagation equation" Eq. (19) becomes 



D 


'( a .\ 2 D <f>' 


4 

n 


'l«t\ iil 


dt 


\ a / [> r 


a 0 


.U/ 



where 

D -3 + A l 

DT ST _ a-. 



( 22 ) 



(23) 



Elrod [1] considered sinusoidal disturbances and found it bene- 
ficial to perform a frequency analysis by representing 4> in the form: 



4 - F(^)e iWY 



(24) 



W is a dimensionless frequency number. 

It was also assumed that 

= l ft) (25) 

The resultant differential equation for n(£) is 

V?) = ° ' ( 26 ) 
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where 



c< =: -2z- 
O-o 





and 



(27) 



(28) 





5^-3 Zzl 3?T-5- 

2 2 2 . 



(29) 



Using the steady state, one-dimensional, isentropic relationships 
for area and temperature as functions of Mach number (see Shapiro [5]), 
the following relationships hold 




' Chiang [2] arrived at the general solution for Eq. (18) by using 
the WKB method. (See Morse and Feshback [6]). 
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( 32 ) 



and are arbitrary constants. Eq. (24) now is written as 

<f>=(-jY e^p[± x- wJ'MdfJ -f c,eap[Aw(r-f)J 

xe*fjj +C 2 exf[-w(r+f)] 



(33) 



X &xf 



aW 



f 



Chiang [2] states that the C-j term with the exp £iw(x-£)~| indi- 
cates an outgoing wave (in the positive r-direction) while the term 
with exp ^iw(x+£)J indicates an incoming wave; the term exp Qiw 
represents a mean convection; and the factors exp f>s/ 9(€)d(j 
represent attenuation of waves. For the case of an infinite conical 
nozzle, with pressure fluctuations impressed at the throat, C 2 = 0. 

The WKB solution is asymptotically valid in the limit of small Mach 
number M and high frequency number W with an error of the order (M /w ). 

Chiang [2] extended the solution to larger Mach number and smaller 

/ 

frequency number W. A summary of this extension follows: 
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Assume 




c , e '^ Z . d ) + c 4 e ‘ W 5 2 ,(!) 



LvJ'f 



(34) 



After substituting into the previous differential equations for n 

and making comparisons, asymptotic approximations of Z-|(£ 0 ) and 

Z 0 (£ ) are determined. 

Z v o 



Z ,(0 = 






fw /%«') if 

oO 



(35) 



£ 1 <*t 

w co 

2,(9 = [-4 £%<«■) J r] 
* TVs') 



(36) 



where £ Q is sufficiently large (i.e., far enough away from the throat) 
so that the WKB method holds. Differentiating Eqs. (35) and (36) 
yields 




(37) 



(38) 



It is further shown that Z-j and are complex conjugates and that 
Eqs. (35) - (38) can be used for initial conditions to solve the 
differential equations for Z-j and Z 
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Let 



z, = u ,+iv; 

Z = U z +Xv L 



then 



U| = u z 






After changing variables by using 



and 



x= M 



* d% 



Chiang showed that 



$ = - 4 -^ 2 . 



and arrived at the following differential equations 



U,"(x) + -f- 



4M- + 2 w ur'(V) 

C? A 



- 4 M. u,(xl = o 



Vr"(*) +-f- v'(x) - 2 W u; (x)l 



~-§M ir-fa) = O 
£4 ' 



(39) 



(40) 



(41) 

(42) 



(43) 



(44) 



(45) 
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with the initial conditions 



u.OO =/ 




u; (^ = 0 




vr ( ^0) ~ v 




vr (^0) ~ 2.w 


fW 



( 46 ) 



where x Q is the value of x corresponding to £ = .3, S, <t>, and g 

have all been previously defined as functions of M and are therefore 
functions of x. 

It is always possible to select a Mach number M small enough to 
satisfy the restriction on the WKB solution, 



<«J 

w 2 - 



(47) 



The amplitude and phase relations between pressure and mass 
flux perturbations are determined in the following manner. Eqs. (8) 
and (10) are used to obtain 

± %L = z ( 48 ) 

9* 2>r 



J_ 24 I - 3±L' + <7. JJil ' (49) 

f 7>* ~ die ^ 
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The mass flux perturbation is written 



H)' = e * %e' 



After introducing the dimensionless pressure and mass flux perturba- 
tions: 





u, 


(51) 


e % 




(52) 



Chiang [2] showed that 



3-r _ V M iz — ? ? 

-2^. m 

JT r 3Y 



(53) 



Furthermore, since both tt and y are periodic in time and contain the 



factor exp £ iwxj 



IT _ 
/* 



3TT 

^-r' 

3r 



(54) 



For the case of an infinite conical nozzle with sinusoidal 
pressure fluctuations impressed at the throat, C 2 in Eq. (34) be- 
comes zero and 
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(55) 



\(tj = C, ©^[-xWf] (u.+iv;) 



therefore 



x e<cp[Aw(r- j)T ( U| +iv;) 

Taking the logarithm of $ and differentiating gives 



(56) 




Substitution of Eq. (57) and C58) Into Eq. (531, together with. 
Eq. (54) yields 



rr . 

a 

where 

N 






~~ 1-^2- 

37^ 




= amplitude factor = 



N, z + V k 

N , 1 + < J 



y ■ phase angle = arctan f ^LJ^i ^4 \ 

V IV, W 3 + N t A/*y 



(.59) 



(60) 



(61) 
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and 






"4" p* /^u^vr 1 ! *cU 



4 . y- dv7_ ) 
<** / 



N t ^ - *(*■* *') ( Ul # ~ ^ 4 t) 



( 62 ) 



N^~N, 



)V^ = - Wa. ~ W -X 



_]_ 

2 . 



A dimensionless parameter designated as the modified Strouhal 
number S* was defined by Chiang [2] 



C* - 63 r 
* ” % 



(63) 



and shown to be related to the frequency number 



w = M s 



p t 



(64) 



For a given Mach number M, it is more convenient to express the 
mass flux-pressure relationships in terms of the modified Strouhal 
number S* rather than the frequency number W. 
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III. NUMERICAL SOLUTION 



Equations (44) and (45) are transformed to a system of four first 
order differential equations by defining the variables: 



y-j (x) = u-j (x) 
y 2 (x) “ (x) 

y 3 (x) = v^x) 
y 4 (x) = v'(x) 

The resultant system of differential equations is 



( 65 ) 



/ 

yi - 



*2- 



^3 



- ~—[d* V* + £ W + * 



5 

f* 



( 66 ) 



yi - 






A numerical solution for this system of equations was obtained 
using Hamming's Modified Predictor-Corrector Method of numerical inte- 
gration. (See Ralston and Wilfe [7]). It is a stable, fourth order 
integration technique which uses four preceding points for the evalu- 
ation of each successive point. PROGRAM TWO makes use of an IBM [8] 
subroutine, DHPCG, to perform the numerical integration. 

Evaluation of the initial conditions v-j (x q ) and v-j (x Q ) is accom- • 

plished using PROGRAM ONE. Subroutine GOPHI evaluates the function 

% 

g(x)/4>. Numerical integration of v-|(x Q ) is performed using a twelve 
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point Gaussian Quadrature formula which integrates polynomials up to 
degree 23 exactly. Subroutine DQG12 is used for this purpose. (See 
IBM [8], p. 299.) 

With u-j , u-j , , and known, subroutine OUTP of PROGRAM TWO 

computes values of amplitude Factor N and phase angle ip for a given 
Mach number M and frequency number W. The corresponding modified 
Strouhal number S* is also computed. 

The programs have been written for the case of air with y = 1.4. 

IV. EXPERIMENTATION 



A. EXPERIMENTAL APPARATUS 

To approximate the characteristics of an infinite conical nozzle, 
a nozzle 17.5 inches in length with an 8 inch diameter entrance and a 
1 inch diameter exit was designed and constructed. It subtended a 
half-cone angle of 11.3 degrees and featured an entrance-exit planar 
area ratio of 64. This ensured that stagnation conditions existed at 
the entrance for all subsonic flows in the nozzle. (See Shaprio [5], 
p. 614.) 

Figures 1 and 2 show the overall arrangement of the experimental 
apparatus and associated electronic equipment. Figure 3 schematically 
depicts the test section. 

Compressed air from two supply tanks passed through a pressure 
regulating valve and into the 10 inch diameter stagnation chamber 
through a baffling arrangement. The stagnation chamber contained 
honeycomb material with a length/diameter ratio of 10 to dampen out 
any disturbances in the plenum air. A fine mesh screen was also installed 
to hasten turbulence decay prior to the entrance of the flow into the nozzle 
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Sinusoidal pressure disturbances in the nozzle exit plane were 
created by a vibrating piston located in the exhaust chamber. The 
piston was driven by a MB Electronics Model EA 1250 Vibramatic Exciter 
powered by a Model 2120 Amplifier using a sinusoidal signal from a 
Hewlett-Packard Model 200CD Wide Range Oscillator. The vibrator was 
mounted on three horizontal tracks so that the zero position of the 
piston could be adjusted for optimum operation. 

Dynamic pressure measurements in the exit plane were recorded with 
a Bruel and Kjaer Model 4138 1/8 inch Condenser Microphone with Adaptor 
UA0036 and Model 2618 1/4 inch Preamplifier connected to a Model 2803 
Two-Channel Power Supply. 

Mass flux (pq) measurements were performed using a Thermo -Systems, 
Incorporated Model 1050 Constant Temperature Anemometer with a DISA • 
Type 55F33 Hot Wire Probe. The hot wire anemometer measured both steady 
state (DC) and dynamic (AC) components of mass flux. The anemometer DC 
voltage was measured on a Simpson Model 2700 Digital Voltmeter. 

The anemometer and microphone signals were passed through identical 
Krohn-Hite Model 3750 Filters (band-pass) to preserve phase. The AC 
voltages were measured by two Hewlett-Packard Model 3400A RMS Voltmeters. 
The filtered dynamic signals were displayed on a Hewlett-Packard Model 
132A Dual Beam Oscilloscope for visual verification of system response. 

Phase angles between pressure and mass-flux fluctuations were meas- 
ured on a Dranetz Model 305-PA-3002 Digital Phase Meter using the output 
signals from the oscilloscope amplifiers. Figure 4 shows a block diagram 
of the electronic measurement equipment. 

Manometers were used to measure stagnation and exit static pressures. 
Stagnation temperature was measured by a mercury-in-glass thermometer. 
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B. PROCEDURE 

Mass flux (pq) calibration of the hot wire anemometer was performed 
using a Thermo -Systems, Incorporated Model 1125 Calibrator. The method 
of least squares was used to fit an eighth order polynomial in the 
quantity (V - V Q ) to the calibration points. (See Computer PROGRAM 
THREE.) Microphone calibration was performed with a Bruel and Kjaer 
Piston Phone: Type 4220. 

Each experimental run consisted of first setting a steady state 
flow rate by maintaining a constant pressure in the stagnation chamber. 
The frequency and amplitude of the pressure disturbances were set on 
the oscillator. Both band-pass filters were set to pass only a narrow 
band at the disturbance frequency with gain settings of 20 db. Data 
was taken over the range of frequencies from 20 Hz to 250 Hz. 

C. DATA REDUCTION 

Reduction of experimental data was accomplished on a digital com- 
puter through the use of PROGRAM FOUR. 

Steady state mass flux was computed by using the anemometer DC 
voltage in the evaluation of the hot wire polynomial. The magnitude of 
mass flux perturbation was obtained by first correcting the anemometer 
RMS voltage and then adding it to the DC voltage to yield a peak voltage. 
The hot wire polynomial was then evaluated for the peak voltage, giving 
a maximum value of mass flux. The steady state value was then subtracted 
from the maximum value to arrive at the perturbation. 

Microphone RMS voltage was similarly corrected to obtain a peak value 
and then multiplied by the microphone sensitivity to yield the magnitude 
of the pressure perturbation. A 6 db. narrow band insertion loss in the 
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band-pass filters was also accounted for in the computations of dynamic 
quantities. 

Mach number was computed from the one dimensional, isentropic 
pressure relationship. The velocity q used in the computation of 
Modified Strouhal number S* was determined by dividing the mass flux 
by a value of p arrived at by using the perfect gas law and the isen- 
tropic pressure-density relationship. 

V. DISCUSSION OF RESULTS AND CONCLUSIONS 



Figures 5 through 14 show a comparison of the experimentally deter- 
mined amplitude factors and phase angles with values predicted by the 
Modified WKB solution for Mach numbers: .084, .117, ,201, .278, and .340. 
Both amplitude factor and phase angle are plotted against modified 
Strouhal number S*. The smooth curves are the modified WKB solutions 
and the experimental results are represented by crosses. 

For low frequency sinusoidal pressure disturbances with correspond- 
ing small values of S*, the phase angle is -180 degrees in accordance 
with the modified WKB solution which asymptotically approaches that 
value. An increase in pressure results in a decrease in mass flux and 
physically it is easy to visualize an effective phase angle between the 
rising pressure and the falling mass flux. Mathematically, this effec- 
tive phase angle is simply 180° plus the value of the actual phase angle. 
Thus, for low frequencies there is no effective lag between u and y. 

The WKB values of amplitude factor monotonical ly increase from unity 
for low S* to the asymptotic limit for extremely high S*. This means 
that the’magnitude of the mass flux response diminishes as the frequency 
of the pressure disturbances increases. 
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Clviang [2] noted that the quasi-steady relationship was probably 
valid for low values of S* but became increasingly inaccurate for larger 
values of S*. 

The experimental values of amplitude factor agree quite well with 
those predicted for all Mach numbers tested. The experimental values 
for phase angle differed from those predicted but exhibited a similar 
trend. For the low frequency tested (20 Hz) all experimental phase 
angles were close to -180 degrees while the WKB values were slightly 
higher. For higher frequencies, the experimental phase angles became 
larger than the WKB predicted values. This apparent phase shift can 
be attributed to effects of viscosity in the flow. Viscosity serves 
to dampen the flow system with a resultant increase in effective phase 
angle over the undamped system, i.e., that assumed in the WKB analysis. 

Physical limitations on the test apparatus prevented investigations 
in the regions above M = .34 and disturbance frequencies of 250 Hz. 
Electronic noise and vibrations generated by the air compressor prevented 
its being run during experiments. Consequently, the air supply in the 
tanks could not sustain runs over a frequency range at higher Mach num- 
bers. The other problem encountered was that of turbulence. Early in- 
vestigations revealed that turbulence in the nozzle was prevalent at 
frequencies below 200 Hz. Inclusion of the fine mesh screen in the 
stagnation chamber helped greatly to dissipate turbulence or at least to 
break it up into high frequency turbulence. Power density spectrums of 
the hot wire anemometer signals, obtained by taking the Fourier Transform 
of the autocorrelation function, indicated the presence of turbulence 
with frequencies above 250 Hz. This high frequency turbulence created 
too high a signal to noise level for meaningful data to be obtained 
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above 250 Hz. A PAR Model 1 01 A Correlation Function Computer and a 
Model 102 Fourier Analyzer were used to determine the power density 
spectrums . 

It is concluded that the modified WKB solution, with its ideal 
flow assumptions, does offer a reasonable approximation of the ampli- 
tude factor for pulsating subsonic flow of air in a conical nozzle, 
as long as the flow originates at near-stagnation conditions. Pre- 
diction of phase angle is conservative and could be used as a lower 
bound . 

It is recommended that further investigations be pursued at higher 
Mach numbers. This would require both a quiet compressed air supply 
and a higher power vibrational exciter. It is felt that the present 
location of the air compressor within the laboratory is unsatisfactory 
and that the compressor should be isolated from the laboratory by 
either moving it outside of the building or by some other suitable 
means . 

It is also recommended that an analysis be undertaken which would 
account for the effects of viscosity upon the dynamics of a pulsating 
flow field in a convergent nozzle. 







Figure 1 . TEST APPARATUS 













Figure 2. ELECTRONIC EQUIPMENT 
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MICROPHONE THERMOMETER 
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FIGURE 3. SCHEMATIC OF TEST APPARATUS 




DIGITAL 

VOLTMETER 




O 2: 
m 



31 



FIGURE 4. BLOCK DIAGRAM OF ELECTRONIC MEASUREMENT EQUIPMENT 





FIGURE 5. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR WITH 
• MODIFIED WKB SOLUTION, M = .084 
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FIGURE 6. COMPARISON OF EXPERIMENTAL PHASE ANGLE WITH 
MODIFIED WKB SOLUTION, M = .084 
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FIGURE 7. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .117 
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FIGURE 8. COMPARISON OF EXPERIMENTAL PHASE ANGLE - 
WITH MODIFIED WKB SOLUTION, M = .117 
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AMPLITUDE FACTOR N 



FIGURE 9. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .201 
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PHASE ANGLE Y 



FIGURE 10. COMPARISON OF EXPERIMENTAL PHASE ANGLE 
WITH MODIFIED WKB SOLUTION, M = .201 " 
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AMPLITUDE FACTOR N 



FIGURE 11. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .278 
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PHASE ANGLE ¥ 



FIGURE 12. COMPARISION OF EXPERIMENTAL PHASE ANGLE 
WITH MODIFIED WK8 SOLUTION, M = .278 
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AMPLITUDE FACTOR N 

FIGURE 13. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .340 
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PHASE ANGLE Y 

FIGURE 14. COMPARISON' OF EXPERIMENTAL PHASE ANGLE ' 
WITH MODIFIED WKB SOLUTION, M = .340 
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